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ABSTRACT 

We investigate the nonlinear ac transport through a quan- 
tum wire with an impurity in the presence of finite range 
electron-electron interactions. We discuss the influence 
of the spatial shape of the ac electric field onto transport 
properties of the system and find that the scaling behavior 
of the occupation probability of the sidebands depends on 
the range of the voltage drop. 

1. INTRODUCTION 

Time dependent quantum transport has attracted a lot of 
interest since the works of Tien and Gordon 111] and Tu- 
cker more recently, theoretical findings and ex- 
periments on quantum dots |^] and on superlattices [Q] 
renewed the interest in photon-assisted transport in semi- 
conductor nanostructures. In particular, the possibility to 
investigate experimentally time-dependent transport th- 
rough mesoscopic systems has opened the way to a deeper 
understanding of new effects strongly relying on the spa- 
tiotemporal coherence of electronic states. Moreover, in 
most time-dependent experiments like electron pumps 
[% U, photon-assisted-tunneling |Q ||, [hJ, and lasers 
[ lip require an analysis going beyond the linear response 
theory in the external frequency. Thus, many efforts have 
been devoted, in last years, to the theoretical investigation 
of nonlinearities in semiconductor nanostructures [|l2[ |l3|] , 
electronic correlations Jl4| , OJ, and screening of ac fields 

The Tien-Gordon formula, according to which the dc 
component of the photo-induced current is given by a 
superposition of static currents Iq (the currents without 
the ac field) weighted by integer order Bessel functions, 
is represented by the following formula 
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shown that the side-band peaks depend on the screening 
properties of the system [|T7h; moreover theoretical inves- 
tigations for superlattice microstructures showed an £1~ 2 
dependence of the transmission probability spectrum of 
the photonic sidebands (that is the argument of the Bessel 
functions), when a nonlocalized (a finite range) ac driving 
was taken into account [ 18 jl], 

In this paper, we investigate how Id electron-electron in- 
teraction, in the framework of the Luttinger model 
[22[ |23| ], nonlinearities, due to the presence of an impu- 
rity, and a finite range ac electric field affect the photo- 
induced current. We will show that the TG formula is 
still valid, but the argument of the Bessel functions is not 
anymore linearly dependent on 

In the time dependent regime the nonlinearity of the sys- 
tem gives rise to frequency mixing and harmonic gen- 
eration. Earlier treatments of the ac transport consid- 
ered voltages, dropping only at the position of the barrier 
p4] , |25| ], and zero range interactions between the elec- 
trons. Here, both of these are generalized to the more 
realistic situation of finite range of both, the electron- 
electron interaction and the electric field. As a matter 



of fact previous calculations [ 26 1 showed clearly that the 



spatial shape of the electric field does influence ac trans- 
port. 

2. MODEL 

The Hamiltonian for a Luttinger liquid of length L (— > 
oo) with an impurity and subject to a time-dependent 
electric field is H = Hq + Hi mp + iJ ac , where 
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The dispersion relation of the collective excitations, 
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the argument of the Bessel functions is linearly depen- 
dent on the ac voltage intensity V\ and on the inverse of 
the driving frequency (or subharmonic) f2. A selfconsis- 
tent theory, based on the scattering matrix approach, has 



depends on the Fourier transform of the finite range inter- 
action potential [ITS]. We assume a 3d screened Coulomb 
potential of range a -1 projected onto a quantum wire of 
diameter d as a^ 1 . The interaction decays exponentially 



and one gets V cc (x) = (Vha/2)e~ a ^ x K with interaction 
strength Vl For a — > oo, one obtains a zero-range 
interaction. 

The tunneling barrier of height Ui mp is localized at x = 



ffimp = CAmp COS (2yfif0( X = 0)) , (3) 

with the phase variable of the Luttinger model 
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The coupling to the external driving voltage yields 

POO 

H^c = e / dxg(x)V(x, t). 



The electric field is related to the voltage drop by differ- 
entiation, E(x, t) — —d x V(x, t), and the charge density 
is g(x) = kp/n + d x 'd(x) j ' \pk ' . The space-time depen- 
dent electric field, E(x, t) = E& c (x) + E a (x) cos (fit), 
such that E a (x) = E\e~\ x ^ a , gives a voltage drop V\ 
| dxE a (x) = 2E\a. The spatial dependence of the 
dc part of the electric field does not need to be specified, 
as only the overall voltage drop, Vq = f_ dxEi c (x), is 
of importance in dc transport | 



3. METHODS AND RESULTS 

The current at the barrier is given by the expectation value 
I(x = 0, t) = (j(x — 0,t)), where the current opera- 
tor is defined via the continuity equation, d x j(x,t) = 
—edtp(x,t). For a high barrier, the tunneling contribu- 
tion to the current can be expressed in terms of forward 
and backward scattering rates which are proportional to 
the tunneling probability A 2 . The latter may be obtained 
in terms of the barrier height Ut by using the instanton 
approximation [|o|. The result can be written in terms of 
the one-electron propagator S + iR [53], 

/>oo 

I(x = 0, t) = eA 2 / dr c~ 5(t) sin R(t) 
Jo 



dt'V eS (t') 
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with 



S(t) + iR(r) = — / — He{a- 1 (x = 0,oj)} 



(1 — cosojt) coth — — h i s'mujT 



where j3 — l/k^T, w max the usual frequency cutoff that 
corresponds roughly to the Fermi energy [f3l|], and the ac 
conductivity of the system without impurity is [E6| 
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Furthermore, the effective driving voltage is related to the 
electric field by [M 



V eS (t) 



dx 



dt'E(x,t')r(x,t - t') 



V Q + — 1*1 cos (ft*- <p z ), (6) 
e 



where r(x, ui) = cr(x, ui)/a(x, oj), \z\ and ip z are, respec- 
tively, modulus and argument of 
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— J dxE a (x)r(x,Q). (7) 



With the above assumptions about the shapes of the driv- 
ing field and the interaction potential one obtains 
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where fc(f2) is the inverse of the dispersion relation and 



A 2 (it, v, w) = 
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In the following, we concentrate on the results for the dc 
component of the current which does not depend on x and 
is directly given by the current at the barrier, for which we 
only need to know only \z\, 



Idc= J2 Jn(M)Io[Vo + n— 



(10) 



The important point here is that the driven dc current is 
completely given in terms of Iq(Vo), the nonlinear dc 
current-voltage characteristic of the tunnel barrier, 

I (V ) = eA 2 ^°dre- s(r) sin R(t) sin (^^) ■ (11) 

Eqs. (|To|), (|llb generalize results which have been ob- 
tained earlier [|lj but without interaction between the tun- 
neling objects, and also for the Luttinger model with a 
zero-range interaction, together with a ^-function like dri- 
ving electric field |Q] . 

For Vq much smaller than some cutoff-voltage V c whi- 
ch is related to the inverse of the interaction range, Iq oc 
V* /9 ~\ This recovers the result obtained earlier for 5- 



function interaction and zero-range bias electric field | 
When Vq V c , the current becomes linear [[32]]. For in- 
termediate values of Vq, Iq exhibits a cross-over between 
the asymptotic regimes with a point of inflection near V c . 
For zero-range interaction, Iq oc V 1 for any Vq. 
Figure |l| shows the currents Iq, I^c and the differential 
conductance didc/dVb as functions of eVo/hfl for g — 
0.9 and g = 0.5 for zero-range of the driving electric 
field. For g = 0.9 one observes sharp minima in the 
differential conductance at integer multiples of the driv- 
ing frequency in certain regions of the driving voltage V\ . 
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Figure 1: Currents Jo, idc and differential conductance 
didc/dVb at zero temperature as a function of the ratio 
eVo/hil for g — 0.9 (top), g — 0.5 (bottom) for values 
ft = v^a, a = 0, and eVi/hvFCt = I (I — 5 dotted, 
£ = 6 dashed, £ = 7 dash-dotted lines). Currents in units 
of HvFa/eRu differential conductance in units of -RjT , 
tunneling resistance i? t = 2?iw 2 nax /7re 2 A 2 . 



These can be understood as follows. When the strength of 
the interaction is not too large, the region where didc/dVb 
is much smaller than 1 is small compared with hfl thus 

for eV « m, dI dc /dV oc (2/g - 1) \eV„ - M| 2 /9- 2 . 
Then, Eq. ( flO| ) yields near eVb = m?ir2 



d/dc 
dV 



1 — J m (|z|) + const 
x \eV - mhn\ 2/9 ~ 2 



J 2 m (\A) 



(12) 



For g > 2/3, this yields for integer m the cusp-like struc- 
tures observed in Fig. For g < 2/3, no cusps oc- 
cur anymore. In addition, the current I^c is depleted so 
strongly and over such a large region of the bias volt- 
ages that the regime of almost vanishing didc/dVb be- 
comes larger than TiQ, and in general no minima near in- 
teger multiples of the frequency exist. As can be seen in 



the figure, the depths of the cusps depend on the driving 
voltage Vi (oc |z|) which can also be understood from of 
Eq. ([j~2|) which shows that the values of the differential 
conductances at the voltages eVb = mhfl are approxi- 
mately 1 - Jl(\z\). 

It is therefore instructive to look into the behavior of | z \ 
as a function of the frequency. Figure ^ shows the scaling 
exponent v determined from 



v = —vfcv 



d log \z\ 
d log ft 



(13) 
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Figure 2: Scaling exponent v of the argument \z\ of the 
Bessel functions as a function of ft, for ranges of the driv- 
ing field (curves from right to left) aa = 10~ 3 , 10~ 2 , 
10 -1 , 1, 10, 10 2 , 10 3 , and g = 0.5. 

We observe a non-universal cross-over between \z\ oc 
ft -1 , the case discussed by Tien and Gordon [|J which 
corresponds to a driving field of zero-range (a — > 0), 
and \z\ oc ft -2 which is obtained for a homogeneous ex- 
ternal field (a — > oo) [Q. Although the behavior of z 
depends strongly on the parameters of the model in the 
cross-over regime, this does not influence qualitatively 
the occurrence of the cusps. Their existence depends cru- 
cially on the finite range of the interaction, and the condi- 
tion g > 2/3. However, by varying \z\, the depths of the 
minima are changed due to the variation of J^d^l). 
Finally, we have demonstrated that the result which has 
been obtained by Tien and Gordon for tunneling of non- 
interacting quantum objects in ID driven by a mono-chro- 
matic field localized at the tunnel barrier remains valid 
even in the presence of interactions of arbitrary range and 
shape, and for an arbitrary shape of the mono-chromatic 
driving field. The central point is that the frequency driven 
current is completely given by a linear superposition of 
the current-voltage characteristics at integer multiples of 
the driving frequency, weighted by Bessel functions. 
The argument of the latter contains the amplitude of the 
driving voltage only linearly but the dependence of the 



argument on the frequency and the range of the driv- 
ing field is determined by its spatial shape. However, 
one can easily identify regions where the dependence on 
the frequency becomes very simple. For a driving field 
which is localized near the tunnel barrier, the integral in 
Eq. ([7]) can be evaluated approximately by noting that 
r(x,SY) varies only slowly with x and can be taken out 
of the integral. Then, \z\ = eV\/Ml which corresponds 
to the result of Tien and Gordon [Q]. In the other limit 
of an almost homogeneous electric field, E\ = V\/a, 
one needs to calculate the spatial average of r{x, f2) [poj]. 
This gives a(k = 0,ft)/a(x = 0,f2) w fi , since 
a{x = 0, Q) « const. This implies \z\ oc J7 -2 . Such 
a frequency dependence has been discussed earlier for 
non-interacting particles [18|. Here, we see that it is 
valid under quite general assumptions also for interact- 
ing particles. A possible method to detect this behavior 
experimentally is to investigate the real part of the first 
harmonic of the current through the tunnel contact and to 
determine the current responsivity which is given by the 
ratio of the expansions of I& c and the first harmonic to 
second and first order in \z\, respectively [g], 
Given the above result for the driven dc-current, the gen- 
eral behavior of the differential conductance as a function 
of eVo/hil can be straightforwardly obtained. Of special 
interest is the occurrence of cusps at eVo /Ml — m(m in- 
teger) which appear to be quite stable against changes in 
the model parameters. A similar result has been discussed 
earlier [p3[], but for a small potential barrier between frac- 
tional quantum Hall edge states which implies zero-range 
interaction. In the case discussed here, the finite range of 
the interaction is crucial for obtaining the cusps, due to 
the absence of a linear contribution towards the current 
for small voltage which is characteristic of tunneling in 
ID dominated by interaction. The cusps could be used to 
frequency-lock the dc part of the driving voltage. 
To summarize, we have shown how the electron correla- 
tion and the spatial distribution of a driving field deter- 
mine the anomalous scaling of the photo-induced current 
and the mode locking patterned structure of the nonlinear 
differential conductance. 
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